use dta\Turnout_District_Level, clear 

foreach year in 1912 1915 1918 1921 1924 1927 { 
twoway (lfitci ch_turnout SMD_margin_pre if year==`year', estopts(cluster(PR_district))) (scatter ch_turnout SMD_margin_pre if year==`year', msymbol(Oh)), ///
title(`year') ytitle("") xtitle("") legend(off) graphregion(fcolor(white)) scheme(s2mono) 
graph save figures\gph\turnout`year'.gph, replace
}

cd figures\gph\
graph combine turnout1912.gph turnout1915.gph turnout1918.gph turnout1921.gph turnout1924.gph turnout1927.gph, ycommon xcommon graphregion(fcolor(white)) scheme(s2mono) 
*graph export ..\turnout_V2.tif, replace
*graph export ..\turnout_V2.eps, replace
graph export ..\Figure8.tif, replace
graph export ..\Figure8.eps, replace
cd ..\..\
